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We formulate a self-consistent non-minimal five-parameter Einstein- Yang-Mills-Higgs (EYMH) 
model and analyse it in terms of effective (associated, color and color-acoustic) metrics. We use 
a formalism of constitutive tensors in order to reformulate master equations for the gauge, scalar 
and gravitational fields and reconstruct in the algebraic manner the so-called associated metrics for 
the Yang-Mills field. Using WKB-approximation we find color metrics for the Yang-Mills field and 
color-acoustic metric for the Higgs field in the framework of five-parameter EYMH model. Based 
on explicit representation of these effective metrics for the EYMH system with uniaxial symmetry, 
we consider cosmological applications for Bianchi-I, FLRW and de Sitter models. We focus on the 
analysis of the obtained expressions for velocities of propagation of longitudinal and transversal color 
and color-acoustic waves in a (quasi)vacuum interacting with curvature; we show that curvature 
coupling results in time variations of these velocities. We show, that the effective metrics can 
be regular or can possess singularities depending on the choice of the parameters of non-minimal 
coupling in the cosmological models under discussion. We consider a physical interpretation of such 
singularities in terms of phase velocities of color and color-acoustic waves, using the terms "wave 
stopping" and "trapped surface". 



I. INTRODUCTION 

The coupled system of Einstein- Yang-Mills-Higgs (EYMH) equations form a mathematical basis for the well-known 
self-consistent model of interaction of gravitational, gauge and scalar fields (see, e.g., [l| for review and basic refer- 
ences). The Einstein- Yang-Mills-Higgs theory unifies two important trends in modern theory of gravity. The first one 
is represented by the so-called Einstein- Yang-Mills model, which can be considered as a non-Abelian generalization 
of the Einstein-Maxwell model. The second trend is connected with the investigations of interaction of gravitational 
and scalar fields. SU(n) symmetric EYMH theory synthesizes the ideas and methods elaborated in both models, 
inherits well-known results and presents some qualitatively new ones. One of the most interesting results obtained 
in the framework of the Einstein- Yang-Mills- Higgs model is connected with the search for exact solutions describing 
non-Abelian field configurations of the monopole type (see, e.g. , B|-[§|). Other important applications of the EYMH 
theory are the modeling of cosmological evolution (e.g., [lfj, [U|) and modeling of particle dynamics in the field of 
non-Abelian monopoles (e.g., [12] )■ Why is the EYMH model interesting for cosmological applications? As it was 
stressed in [10( , numerous models of inflation in the Early Universe involve into consideration the multiplet of Higgs 
scalar fields <fr. Higgs scalar fields are coupled with the gauge Yang-Mills field, characterized by potential four- vectors 
Afc, the so-called minimal coupling being realized by means of the gauge-covariant derivative D m & = <9 m «& + [A.,„, «&]. 
Since contributions of the Higgs fields into the total stress-energy tensor of the system are considered in many models 
as the dominating ones at the inflation stage, the gauge counterpart of the Higgs field, the Yang-Mills field, should also 
be included into the master equations. As for t he p resent stage of the Universe evolution, the models with Higgs fields 
appear, first, in the theory of dark matter [l3l Il4j. second, in the theories of dark energy [l5|, [l6|. Just the discovery 
of accelerated expansion of the present Universe revived an interest in non-minimal Einstein- Yang-Mills-Higgs models 
in the context of search for an explanation of the dark energy phenomenon. One of the most advanced trends of 
such type is represented by the so-called f(R) theories (e.g., [17|). which are based on the appropriate (nonlinear) 
modifications of the Einstein-Hilbert part of the total Lagrangian (R/8ttG) linear in the Ricci scalar. Another version 
of non-minimal extension of the EYMH theory is characterized by the non-minimal modification of a scalar field 
contribution to the total Lagrangian (e.g., [Hj]). Our approach to the non- minimal modification of the EYMH theory 
is presented below. 

The discussion concerning the non-minimal coupling of gravity with fields and media started at the end of the 
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60th. The historical details, review and references related to the non-minimal interaction of gravity with scalar and 
electromagnetic fields can be found, e.g., in [2(J HH . As for the generalization of the concept of curvature coupling 
for the case of gauge field, there are two different ways to establish a non-minimal Einstein- Yang-Mills theory. The 
first way is a direct non-minimal generalization of the Einstein- Yang-Mills (EYM) theory containing derivatives of the 
second order at most [22J . In the framework of this approach Miiller-Hoissen obtained the non- minimal EYM model 
from a dimensional reduction of the Gauss-Bonnet action [23| • We follow the alternative way, which is connected with 
a non-Abelian generalization of the non-minimal Einstein-Maxwell theory along the lines proposed by Drummond and 
Hathrell for the linear electrodynamics [24|. Based on the results of the paper [25} . we considered in [2j]|, |2jJ |28( a three- 
parameter gauge-invariant non-minimal EYM model linear in curvature. Next natural step was a formulation of a 
non-minimal Einstein- Yang-Mills- Higgs (EYMH) theory, and this process, of course, also admits different approaches. 
Taking into account the Higgs scalar multiplet we follow, first, the ideas, proposed in [29|, HO, HH HH, second the 
concept of derivative coupling, developed in [H, Hi], [H[ . 

In this paper we establish a five-parameter non-minimal Einstein- Yang-Mills-Higgs model. The first three coupling 
parameters, qi, q2, 93, describe a non-minimal interaction of Yang-Mills field and gravitational field. The fourth 
and fifth parameters, 54, 55, describe the so-called gauge- invariant non- minimal "derivative coupling" of the Higgs 
field with gravity. Since the gauge-invariant derivative, D m <&, contains the potential of the Yang-Mills field, the 
corresponding non-minimal term is associated with "triple" interaction, namely, gravitational and scalar fields, gauge 
and scalar fields, and gauge and gravitational fields. 

It is clear that curvature induced effects predicted by the non-minimal EYMH theory should be examined by the 
analysis of particle dynamics and wave propagation. Wave propagation in media is known to be described in the 
framework of effective metric formalism (see, e.g., [lHHS HH for the history details and references). The first metric 
from the class of effective ones, namely, the optical metric, was introduced by Gordon [l9j| for isotropic media. In the 
effective spacetime with optical metric light propagates as in vacuum, light rays follow geodesies [40J and the wave 
vector is a null vector of this metric. Written in terms of the optical metric, the constitutive equations, linking the 
excitation tensor and the Maxwell tensor in electrodynamics of isotropic media, have the same formal structure as 
for vacuum electrodynamics [41], [42|] • The second representative of the class of effective metrics is the acoustic metric 



36, 37, 38]. Using this metric, one can regard sound waves as quasi-particles moving in the effective spacetime (e.g., 
II). 



The covariant Fresnel equation is known to be a straightforward way to obtain effective metrics in electrodynamics 
(see [42j and references therein). It appears as a product of the geometrical optics approach and can be reduced to one 
eikonal equation in the isotropic case. Another way is connected with algebraic analysis of the constitutive equations. 
In electrodynamics this approach is developed and successfully used by Hehl and Obukhov (42j. In particular, these 
authors obtained the Minkowski metric as a kind of "tensorial square root" of the fourth-rank constitutive tensor, 
linking the excitation tensor and Maxwell (field strength) tensor in the so-called premetric electrodynamics. 

When the medium is anisotropic, light propagation is accompanied by the phenomenon of birefringence 1441] . When 
anisotropy is of the uniaxial typ e the Fresnel equation introduces two effective (optical) Lorentzian metrics [431 > which 
define two light cones [H, Q E^|. In the biaxial anisotropic media the effective (optical) metrics are in general 
non-Lorentzian [45j , but again the birefringence revives bi-metricity. 

Another way (algebraic) yields the same result. As it was shown in [461 ] . the fourth-rank constitutive tensor can 
be generally reconstructed out of two symmetric second rank tensor fields, thus introducing the so-called associated 
metrics. Besides, it was shown in (46| . that the transition between different representations of the constitutive tensor 
in terms of different pairs of associated metrics is governed by invariance properties in an associated two-dimensional 
internal vector space. Despite the fact, that interpretation of these associated metrics as Lorentzian optical metrics 
is restricted to the uniaxial case, the representation of constitutive equations in terms of associated metrics is valid in 
linear electrodynamics of an arbitrary anisotropic medium (exceptional case relates to the presence of skewons [421] ) . 

Here we suggest the generalization of the formalism of associated metrics for the case of gauge field and introduce 
a new type of effective metrics, the color metrics. In order to motivate this generalization, we consider here a specific 
sort of vacuum field configuration, indicated as vacuum interacting with curvature, or equivalently, non-minimal 
vacuum. The five-parameter non-minimal EYMH model, considered below, is a convenient model for our purpose due 
to two reasons. First, the non-minimal EYMH model admits a clear Lagrangian formulation, the master equations 
for gravitational, gauge and scalar fields being obtained by a direct variational procedure. Second, the non-minimal 
master equations look like the corresponding equations for the anisotropic medium. In other words, in the presence 
of the curvature coupling of the gravitational, gauge and scalar fields the vacuum non-minimal EYMH model can be 
reformulated as a minimal one, but in some effective medium (quasi-vacuum, or vacuum, interacting with spacetime 
curvature). As a result, the effective metrics for the Yang-Mills field (color metrics) and for the Higgs field (color- 
acoustic metric) appear in a natural way. 

The paper is organized as follows. In Section [IT] we introduce basic formalism of the EYMH theory. In Section 
IIIII we obtain the non-minimal extensions of the master equations for the gauge, scalar and gravity fields in the 
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framework of five-parameter model. In Section IIVI we discuss the formalism of multi-metric representation of the 
constitutive equations for the Yang-Mills and Higgs fields, and apply this formalism to the case of uniaxial non- 
minimal vacuum. In Section |V] we study in detail the cosmological applications of the effective metric formalism for 
Bianchi-I, Friedmann-Lemaitre-Robertson- Walker (FLRW) and de Sitter models. Conclusions summarize the results. 
In Appendix [A] we consider the symmetry properties of the non-minimal constitutive tensor for the Yang-Mills and 
Higgs fields. Appendix [B] contains the WKB-analysis of the problem, which explains, why the associated metrics 
for the Yang-Mills field can be treated as color metrics. In Appendix [Cl we present the tidal force, which gives an 
alternative description of the particle motion in the presence of non-minimal (curvature induced) interactions. 

II. PREAMBLE 

Our aim is to consider the EYMH model with the action functional of the following type 

S(eymh) = J d^x^—g |^±^ + \C^F$F%& - Cg )W A* (o > + n* 2 )} • (1) 

Here g = det(gik) is the determinant of a metric tensor g^., R is the Ricci scalar, A is the cosmological constant, 
the multiplet of real tensor fields F^' describes the strength of gauge field, the symbol <E>( a ) denotes the multiplet of 
the Higgs scalar real fields, A denotes the gauge- invariant derivative, V(<& 2 ) is a potential of the Higgs field, and 
<f> 2 = (<j>( a )<j>( a )). Latin indices without parentheses run from to 3, (a) and (6) are the group indices, the summation 
over repeating indices is assumed. The quantities C 1 ^™^ and C(a)(b) denote the so-called constitutive tensors for 
the gauge and scalar fields, respectively. They contain neither Yang-Mills strength tensor F^\ nor gauge-covariant 
derivative of the Higgs field L\<E>( 6 ), i.e., the Lagrangian in (fTJ) is quadratic in these quantities. Since the tensors 
Qa)(&) an< ^ ^(a)(b) con tai n group indices (a) and (b), the metric in the group space G( a )(&), the scalar fields $' a ' and 
some additional directors in the group space, q^ a \ can be used to reconstruct these tensors. The symmetry properties 
of the tensors Ci*V}R and O^um an d their general decompositions are considered in detail in Appendix [A] For the 
U(l) symmetry such theory is discussed in [H|,[47|], here we consider SU(n) (n > 1) models. Let us mention that 
in this paper we restrict ourselves by the condition that only the metric tensor in the group space, G/ )m, is used 
to construct the tensors C^™^ and QK, b y Other constructions, for instance, $(„)$(&) tensor, etc., are considered in 

E3- 

In addition to the real fields F^\ and (the Yang-Mills field potential) we use below the symbols Fife, 
$ and Aj. In literature there are few alternative definitions of these quantities, despite all versions give exactly the 
same result in terms of multiplets of real fields. In order to avoid the ambiguity, we stress, that in this paper we follow 
the definitions of the book [49( (see Section 4.3.) Thus, we consider the Yang-Mills field F m „ and the Higgs field <& 
taking values in the Lie algebra of the gauge group SU(n) (adjoint representation): 

F mn = -igt (a) F^, K m = -iQt {a) At\ * = t (a) $( Q ). (2) 

Here t( a ) are the Hermitian traceless generators of SU(n) group, thus, <fr is considered to be Hermitian, but F mrl and 
Ai are anti-Hcrmitian. The group index (a) runs from 1 to n 2 — 1. The scalar products of the Yang-Mills and Higgs 
fields (indicated by bold letters) are defined in terms of the traces of the corresponding matrices (see, [Hj]), the scalar 
product of the generators t( a ) and t( fc ) is chosen to be equal to: 

!*(&)) = 2Tr t( a )t(6) = G( a )( fc ) . (3) 

The symmetric tensor G( a )( b ) plays a role of a metric in the group space and the generators can be chosen so that the 
metric is equal to the Kronecker delta. The representation @ , (J2J) allows to consider the multiplets of the real fields 
{Fmn}, {A^} and as components of the corresponding vectors in the n 2 — 1 dimensional group space. The 

operations with the group indices (a) are assumed to be the following: the repeating indices denote the convolution, 
and the rule *&( a ) = G '( a )(b)^ f° r the indices lowering takes place. In such terms the gauge invariants in the action 
functional (H2J) is reduced to 

4> 2 = (*, *) <^ ia} ^ , (F mn , F mn ) -G 2 F^F^ , 

(An*, D m *) => b m ^ ia) D m ^ { a) ■ (4) 
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The Yangjylills fields F^l are connected with the potentials of the gauge field by the well-known formulas (see, 
F m „ = V m A„ - V„A m + [A m , A„] F^l = V m A^ - V n A$ + Gf%,,A^A^ . (5) 



Here V m is a covariant spacetime derivative, the symbols f^n^ic) denote the real structure constants of the gauge 

group SU(n). The gauge covariant derivative Z? m <f? = t( a )Z) TO $( a ) is defined according to the formulas (|43]> Eqs.(4.46, 
4.47)) 

D m * = V m * + [A m , *] =► D m $ (o) = V m $ (a) + 0/ ( $ (c) 4?* (c) • (6) 
For the derivative of arbitrary tensor defined in the group space we use the following rule [51( : 

The definition of the commutator in §5§ and ([S]) is based on the relation 

[t( a ),t (h) ] = i/ ( ( c i )(b) t (c) , (8) 

providing the formula 

/(c) (o) (5) = G (c)(d )/ ( ( rf o)(6) = -2i Tr [t( ),t(6)]t (c ). (9) 



The structure constants fr a )(b)(c) are supposed to be antisymmetric under exchange of any two indices |49i . l50l . l5l| . 
Metric G( a )(b) and the structure constants f ?a)(c) are su PP ose d to be constant tensors in the standard and covariant 
manner [5l| . This means that 

dmG(a){b) = ' ^m%)(ii) = , 9m/ (°j,)( c ) = , An/ (&)( c ) = • (10) 

Furthermore, when the basis t( a ) is chosen to provide the relation G( a )(b) = <5( a )(6) ; it holds : 



;/ < (!)( c )/ < (#) - 5 (a)(6) - ( n ) 



n 
and 

(t(o) , t(b) } = t(„)t( 6 ) + t( 6 )t( ) = — <5 (a )(b)I + rf ( (i )(b) t (c ) (12) 
[5l| with the completely symmetric coefficients d( c )( a )(&) (I is the matrix-unity). The tensor FfK satisfies the relation 

' ( „i - " 

the asterisk introduces the dual tensor 



D k *Ff*=0, (13) 



*Tpik ^ Akls rp f-\ A\ 

where e lkls = -^=E lkls is the Levi-Civita tensor, E lkls is the completely antisymmetric symbol with E 0123 = — -E0123 = 
1. 9 

As a first step we consider here the so-called non-minimal EYMH model, for which the constitutive tensors are 
constructed using the spacetime and group metrics, the Riemann and Ricci tensors, and Ricci scalar. Other possibilities 
will be studied in future papers. 
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III. FIVE-PARAMETER NON-MINIMAL EYMH MODEL 
A. Non-minimal action functional 

In the context of this paper we assume, that the constitutive tensors and ClK, b j, included into the action 

functional |T]), contain neither contributions from the gauge field (A^ a \ Fmn) nor contributions from the Higgs field 
($( a ) 5 Vfc$( a ', etc.). The tensors C^™^ and can be generally reconstructed using tensor quantities of several 

types. The quantities of the first type do not contain derivatives of the metric; they include the metric <?j/. itself, 
timelike four-vector U k of the macroscopic velocity of the Yang-Mills- Higgs system as a whole, spacelike four-vectors 
"D k a \ directed along anisotropy axes, tensors in the group space, such as G( )(b), J n,\i c y ■■■ etc. Tensor quantities of 
the second type contain the first derivative of the metric only, they can be rewritten in terms of covariant derivatives 
V m U k and V m T> k a y The tensor quantities of the third type consist of second derivatives of the metric and can be 

represented in terms of the Riemann tensor R l kmn and its convolutions: the Ricci tensor Rf. n = R™ mn and the Ricci 
scalar R = R k . The quantities of the next type may include the covariant derivative of the Riemann tensor VjRikmn, 
the tensor R lkmn 'V m U n , etc. When the terms containing the velocity four-vector XJ % and its derivatives are taken 
into account, one deals with the so-called dynamo-phenomena (see, e.g., [HJ for dynamo-optical effects). When the 
four- vectors T^ k a ^ and its derivatives are inserted into the Lagrangian, one assumes that the system contains a spatially 
anisotropic material subsystem. One refers to the effects induced by curvature, when the Lagrangian contains the 
Riemann tensor and its convolutions, and there are no contributions from U k and T>,, \ . 

In this paper we focus on effects induced by curvature and consider the Lagrangian of the EYMH theory linear in 
the Riemann tensor. The ansatz of linearity with respect to Riemann tensor leads us to the five-parameter model 
with the action functional written in the following (generic) form 

Vmeymh) = J d A x^—g + \F^F^ - D m ^D™<Z> {a) +mH 2 

+ \n lkmn F^F mn(a) - ^t mn D m ^D n ^ a) \ , (15) 

where the so-called susceptibility tensors TZ lkmn and 3? m ™ are defined as follows: 

K lkmn = —R {g lm g kn - g m g km ) + ^{R tm g kn - R m g km + R kn g" n - R km g m ) + q 3 R lkmn , (16) 

5R m ™ = q 4 Rg mn + q 5 R mn . (17) 

This action describes the five-parameter non-minimal EYMH model, and q%, 02, ■ ■ • 05 are the constants of non- minimal 
coupling. The action (fT5")) is the particular case of |T]) with 

G (a){b) , (18) 

and 

Cr a )( b )= G (a)(b)(9 mn + ^ mn ) ■ (19) 

Below we indicate the tensor 

g lk = g ik + W k (20) 

as a color-acoustic metric. 
REMARK 1 

The non- minimal susceptibility tensors n ikmn 02J) and jft m " (fT7)l can be rewritten in terms of irreducible parts of 
the curvature tensor 

n ik m n =Xi gikmn +X2S ikmn +X3C ikmn ^ sjjmn _ ^mn + X 5 S mn (21) 

where C lkmn is the traceless Weyl tensor 

Qikmn _ j^ikmn gikmn Qikmn qtyi q (22) 



r-iikmn 



(A - g m g ) + K 



y ikmn 
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and 

gikmn _ jj_ ^gim gkn gin gkm _^_gkn gim gkniginy^ gmn — ^mn ^LJ^g 



Q lkmn = —R (<f 



-fJ 



(23) 



(we use the notations from the book [531]). The phenomenological parameters Ai, A2, A3, A4 and A 5 are connected 
with qi, q 2 , q 3 , q± and q$ by the linear relations 



Ai = 6qi+3q 2 +q 3 , A 2 = (72+93 , A 3 = q 3 



A 4 = <74+^<75 , A 5 = q 5 . 



(24) 



The representations (TITd]) and (jTTJ) are equivalent to (j2"Tj) from a computational point of view, nevertheless, we use 
below the decomposition (fT6|) taking into account some historical motivation (see, e.g., reviews fl9l. I20I] ) . In particular, 
the phenomenological model proposed by Prasanna in [54| contains the Riemann tensor only (qi = q 2 = q^ = q$ = 
0) and gives a first example of the violation of the strong equivalence principle in the framework of non-minimal 
electrodynamics. The first non-phenomenological model elaborated by Drummond and Hathrell in [24j | is based on 
one-loop corrections to quantum electrodynamics in curved spacetime; the coupling parameters are calculated directly 
in [24[ yielding qi = —5q, q 2 — 13q, q 3 — ~2q (or equivalently Ai = 7q, X 2 = 1 lg, A3 = — 2q). The positive parameter 

<7, q = is constructed by using the fine structure constant a, and the Compton wavelength of the electron A c . 
The models investigated by Horndeski [H, [55[ , by Miiller-Hoissen [23[ , by Miiller-Hoissen and Sippel [56[ , contain 
master equations of the second order only; they are characterized by the following values of the coupling parameters: 
qi = —q, q 2 = 2q, q 3 = —q (or equivalently Ai = — q, X 2 = q, A3 = — q) (here parameter q is arbitrary). These relations 
can be recovered by the ansatz that the non-minimal susceptibility tensor TZikmn is proportional to the double dual 
Riemann tensor *R* kmn , i.e., TZikmn — 7 *R*km-r,i f° r some 7 j2||. Analogously, one can use the Weyl tensor Cikmn 
in the relation TZikmn = ^ Cikmn, for some to [25j |. Thus, there are different physical and mathematical motives to 
specify coupling constants, but here we prefer to consider generic q\ 1 q 2l ... q§, introduced phenomenologically. 



REMARK 2 

The Lagrangian of the non- minimal model linear in curvature can contain also the term £i?<I>( a - )| I>( a ) without gradients 
of the Higgs fields, the corresponding model with £ 7^ is the six-parameter one (see, e.g., [ll[). Here we focus on 
five-parameter EYMH model only. 

REMARK 3 

Let us summarize the assumptions, which we made in the formulation of the five-parameter EYMH model. 

(i) The Lagrangian of the EYMH model is quadratic in F^l and D m ^ {a) (see flU). 

(ii) The Lagrangian of the EYMH model is linear in the curvature tensor (see (|15 )) -(|17 p ). 
(Hi) The group indices in (fl~5|) are coupled by the metric Gt a -\a,\ only. 

The latter assumption leads to an important consequence: it provides all the constitutive tensors to be multiplicative, 
i.e., they are products of the metric in the group space and the corresponding tensors defined in the spacetime. In its 
turn, this consequence guarantees, that for the five-parameter EYMH model the constitutive tensors are symmetric 
with respect to transposition of the group indices. 

B. Non-minimal extension of the Yang-Mills equations 

The variation of the action 5'(nmeymh) with respect to the Yang-Mills potential A^ yields 



D k H\% = -g(Dk<S> {b) )f(a)( b){ c)<S> {c) (g lk + W k ) , 



where the tensor W*\ is defined as 



n_/ik 

«(«) 



__^girn gkn girigkm^ _j_ j^ikmn 



G(a)(b)Fml ■ 



(25) 



(26) 



Equivalently, one can write 



D k H lk = V k H lk + [A kl H ik ] = Q 1 (g lk + 3? lfc ) L> fc #, * 



(27) 
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where 

H lfe = F ik + K tkmn F mn . (28) 

By analogy with electrodynamics of continuous media the tensor H lfc in the non-minimal Yang-Mills equations can 
be called excitation tensor. Below we use the term "color excitation" . Equation (|26| is in fact a linear constitutive 

law, connecting the color excitation tensor H)K and the field strength tensor Fmn (see, e.g., [42l. [lil |57| for the U{\) 
symmetry). It can be clearly rewritten as 

demonstrating that the quantity H lkmn G ( a )(b) pl avs the role of color susceptibility tensor. The term 
g2 fgik _|_ sjjifc^ £>k&, * plays, respectively, the role of color current induced by the Higgs fields. 

C. Non-minimal extension of the Higgs field equations 

The variation of the action iS(nmeymh) with respect to the Higgs scalar field $w yields 



(g mn + $t mn )D n <f> =-m J &. (30) 



By analogy with thermodynamics the four- vector D n G> can be indicated as color four-gradient of the Higgs field (the 
analog of the temperature three-gradient), and the four- vector \E' m can be called the color flux four- vector (the analog 
of the heat flux). Thus, the relation 

* TO = (g mn + ^i mn ) £>„* (31) 
is in fact a non- minimal constitutive law, generalizing the Fourier law in thermodynamics, its general form being 

*{S)=C5ft 6) (A,*) (6) . (32) 



D. Master equations for the gravitational field 

In the non-minimal theory linear in curvature, the equations for the gravity field related to the action functional 
S'(nmeymh) take the form 

-,(NMYMH) 



R*k - ^Rg* = Ag lk + *T£ " 1 a ' ■ (33) 

The principal novelty of these equations in comparison with the well-known equations for non-minimal scalar field is 
associated with the third, fourth, etc., terms in the decomposition 

rp(N MY M H) _ rp{YM) T (H) „(!) -(II) -(III) -(IV) -(V) 



The first term t£ M) : 



rp(YM) — 1 771(a) jpmn T?( a ) -p n ('Qc:^ 

1 ik = jfl^mn^a) ~ *in *k{a) > \ 6b > 

is a stress-energy tensor of pure Yang-Mills field. The second one, T^ k , 

if = £>i$Mt> k $M - \grkDrrM a) D m <S> { a) + 9rk (36) 

is a stress-energy tensor of the Higgs field. The definitions of other five tensors are related to the corresponding 
coupling constants q\, q%, . . . , q§: 

T$ = RT% M) \RikFWF$ + \ [biDk ~ gik&Di] [F&Fffi] , (37) 



8 



T (II) - --a u 

1 ik — 2^ 



A (i 



/) ) „J),(F-WF , „ w U m r»WF 1 n(l) 



pin(a) (^Ji u F kn ^ + R k iF in ( a) ) - R mn F^F kn ( a ) - -D m b m (^F^F k " a] 



-Di 



D, 



(F^F^+b k (f^f^ 



(38) 



rp(IH) 1 nmnls Tj(a) tj ^ Tpls(a) ( jp n p p 



2 F) rri F) r i 



p. n ( a ) p 



k(a) ' r k r i(a) 



(39) 



(g lk b n b n - bibkj [b„M a) b m $> 



'{a) 



(40) 



t' 7 ) — n *(°) 



R?D k <S> [a) + RZDi$ (a) 



1 



(«) 



-^gikD m D n 



(a) 



-D m { D, 



A»* (o) A*(„) 



bM a) b k $ ia) 



(41) 



The Einstein tensor Gik = Rik — ^gikR is the divergence-free one, thus, the tensor p^ MYMH ) j n the right-hand-side 
of (13"3"1) has to satisfy the differential condition 



krp(NMYMH) 



V fe T. 



■ik 



0. 



(42) 



One can prove that it is valid automatically, when F^ is a solution of the Yang- Mills equations (|2T|) , and satisfy 
the Higgs equations (|30| . In order to check this fact directly, one has to use the Bianchi identities and the properties 
of the Ricmann tensor: 



V ' iRklmn ~1~ V ' iRikmn ~1~ V ' kRlimn ) Rklran Rmkln ~\~ Rhnkn 7 

as well as the rules for the commutation of covariant derivatives 

(V;V fe - VkVjM* = A m Ri mlk , 



(43) 



(44) 



(this rule is written here for vectors only). The procedure of checking is analogous to the one, described in [5a | and 
we omit it here. 
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E. Interpretation of the gravity field equations in terms of gravitational excitation and strength 

Since the Riemann tensor R 1 klm can be considered as a gravitational field strength, an analog of Fi k in electrody- 
namics (see, e.g., (59|). one can emphasize an interesting analogy with electrodynamics. Let us mention that one of 
the ingredients of electrodynamics of media is a linear constitutive law 



(45) 



where H lm is an excitation tensor, C lmls is a constitutive tensor and H. lm is a spontaneous polarization-magnetization 
tensor, which does not depend on the Maxwell tensor - The Maxwell tensor and the excitation tensor satisfy 
equations 



V t F kl + V t F lk + V k F u = , V k H lk = --P , 

c 



(46) 



where P is a four- vector of the current of unbounded charges. In analogy with (|45[) we can introduce a linear 
constitutive equation for the gravity field 



7 — i-/ x- r pqls r 

^irakn — r^imkn > ^imkn ^pqls 



(47) 



where Zi mkn is a tensor of gravitational excitation, C imk pqls is a constitutive tensor, and TLi m ki is a tensor of gravi- 
tational polarization. In analogy with ([46]) the field strength R pq i s satisfies the Bianchi identity ((43|) . The Einstein 
equations with look like the convolution of the constitutive equations (|4"T]) with g mn , if we put 



p pqls 1 

imkn 2 



c 



CM [is]) 



/ k 



'mn C-i n 



(MM), 



.,.+/• (MM) £ (NM). 

Jmk 1 >~mn '"mk i> 



_J_ / _ \ rt f (\pq\[lsY) 

a \9ik9mn QinQrak) <7 '-'rt 1 



n pqls _ qs 

'-ik — y 



SfS[ - ~g*> l g ik 



+ ^9 qs { M 



1 



+92 



--au F p{a) F l] + F y (b p F {a) + S p F {a) \ + F p F (a)l 



(48) 



for the constitutive tensor, and 



nj _ (nj{G) nj(G) 1 nj{G) Is 

iLimkn — « I <?fe]m — 'l m [„5fe]i — ^ is 5 5i[nff/c]r 



(49) 



?4 G) = A.9 lfc + r, 



i k 



D i D k -g ik D l D l 



p(a) jpmn 
M mn^ (a) 



2 92 1 9ik An A ( 



pmn{a) pi 



n(a) 



D m D r 



( z?( a ) z? n 
I m fc (a' 



-A 



D 



( p( a ) pin 



p(a) pi 
in (a) 



)]} 
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771 ™0) r m , n 
t k{a)+ t k 



i(a) 



i (a) 



94 



(da 



D n D n - D,D 



D m <f>( a) D m <S> 



(a) 



L> m $ (a) M 



(a) 



(a) 



A^ (a3 -Dfe$(a)]} 



^ik^m^n 



Z) m $ (a) Z)™$ (a) | 



(50) 



for the gravitational excitation. Here we use the following standard definitions for symmetrization and antisym- 
metrization of indices, respectively: 



Finally, the identity 



2 (^mn K^n 



(51) 



-,9ik9 



(NMYMH) 



A<? lfc =0- 



(52) 



written as 



V fc (Z lmkn g mn )=0, 



(53) 



is an analog of the second Maxwell equation with vanishing current. Clearly, the constitutive tensor (|4"5)) entering 

the constitutive law (|47[) . is built using metric, quadratic combinations of the tensor F-£' and color gradient of the 
Higgs field, and does not contain the derivatives of the spacetime metric. 

Thus, in the presented five-parameter non-minimal model three constitutive tensors, C l , (fT5|) , C^. (JTHJ) and 



pqls 



(|48[) , appear in a natural way. It can be considered as a motivation for the construction of more sophisticated 



models with the action functional ([TJ) . 



IV. MULTI-METRIC REPRESENTATION OF THE CONSTITUTIVE TENSORS 

A. General formalism 

Our final purpose is to reconstruct effective metrics (associated, color and color-acoustic) for the five-parameter 
EYMH model with the action functional (fTS"]) for the one-axis Bianchi-I model, FLRW and de Sitter models. In order 
to explain the procedure of their finding for this five-parameter model, we consider, first, a general formalism (see next 
subsection and Apendices), then we reduce basic formulas to the case of uniaxial symmetry, and, finally, we obtain 
desired effective metrics. Let us mention, that we do not introduce new assumptions in addition to the three ones 
fixed above in the Remark 3, but we made a few simplifications, which follow from the symmetry of the corresponding 
cosmological model. For instance, the spacetime symmetry of the one-axis Bianchi-I model requires the symmetry of 
the constitutive tensors to be uni-axial with vanishing cross-effects. This is a reason, why we consider the appropriate 
reduced formulas instead of the general ones. 

1. Decomposition of the tensor Cffiub) 



Let g^™\°^ be tensor fields symmetric with respect to spacetime indices, i and m. Here a — 1,2 is the field 
number, and (a) is the group index (color number), taking only one value for U(l) symmetry and running from 1 to 
n 2 — 1 for the SU{n) symmetry. This Latin index is placed in parentheses. We restrict ourselves by the models with 



C, 



= c, 



ik'ii 



(a)(6) — w (6)(a)' ^ US ' ^ e m ost general representation of this constitutive tensor in terms of 



im(a) . 



IS 



s-iikmn 
°(a)(6) 



1 

2fi 



(c)(d) 

(a)(/3)(a)(b) 



im(a) kn{{3) 

9(c) 9 (d) 



in(a) km{j3) 

- 5(c) 9 (d) 



(54) 



(a)(/3)(c)(d) 
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where the following symmetry of indices is assumed: 

r (c)(d) _ r (c)(d) _ r (d){c) / c c \ 

U ( Q )(/3)(a)(b) - U (a)(/3)(o)(a) ~ U (,3)(a)(a)( ) ' K 00 ) 

the sum being over all possible combinations of field numbers and color numbers. The factor /t is introduced, as in [46| . 
for convenience. The C 1 ^™" tensor has the general form (|A4|) (see Appendix |A"|) . It contains 21 ^" ~% components. 
There are 10 components of gf^ for each a and (a). Thus, when we deal with two metrics, i.e., a = 1, 2 for n > 1 
one has 20(n 2 — 1) field components flS , as well as |(n 2 — \) 2 n 2 {2n 2 — 1) parameters (functions) G^l^/w^ , as a 

tool for reconstruction of 21 ^" components of the CV^^ tensor. The reconstruction of C^ k "^ along the line of 

(|54p is always possible since the corresponding algebraic system is underdetermined. Let us mention that U(l) case is 
a special one, and (see, e.g., [46j) one has 20 field components g lm ( a \ as well as 3 parameters Gi a \(m to reconstruct 21 

components of C lkmn . In the case of SU (2) symmetry we have, generally, 60 field components , 126 components 

of G^\?p\/ w&) and 126 components of C^"™^. Thus, 60 parameters are arbitrary. 

The quantity ^(a)^g)( a )(M i s considered to be a twice covariant and twice contravariant tensor in the n 2 — 1 dimen- 
sional group space, the tensors g7*\ are the vectors in that space. When a, i and m are fixed, the standard unitary 
transformations are assumed to conserve the scalar product G( a )( b ) gj^ a ^ ■ As well, when the spacetime indices 

and group indices are fixed, there exist a linear transformation from one set of g™ l ^ a ^ to another one, and g7^ can 
be regarded as vectors in some effective two-dimensional space [46[ . When we deal with the so-called "parallel fields" 
[ooj ] . i.e., when 

5^ (Q) = ?(a)5 im(Q) , G W(6) *(.)9(6) = 1, (56) 

we can put 

G (aK«(a)(6) = G { a )(l3)G(a)(b)G {c)id) , (57) 

providing C^™^ to have a multiplicative structure, C?^w5 = C lkmn G^Q,y Thus, the model of parallel fields 

for the 5J7 (n) case covers the E7(l) model considered in [40|. When, on the contrary, we assume C 1 ^™^ = C lfcm "G( a )(f,), 

the tensor Gj^Sjum is not obligatory multiplicative, nevertheless, the parallel fields give one of the solutions of the 
reconstruction problem. 

As any symmetric tensor, the quantities can be decomposed with respect to their components parallel and 

orthogonal to the four- velocity U l , 

g ^ = B[^lPU k + V l ffu k + Vffiu* + Sl k a \ a) , (58) 

where 

®{a) - 9(a) U i U k, ^ (a ) = A t ff( a ) U k, = A l 5 (a ) A fc • l° y i 

In the application of this formalism to cosmology (see below) , one obtains directly that 

X>gf = =* 5 * W =C/ i C7 fc +5g Q) =► U i9 ^ a) = U m , (60) 
which means that the velocity four-vector U k is the eigenvector for all tensors fl£*\ with the eigenvalue equal to one. 



2. Decomposition of the tensor C(*)(m 
We assume that Qa)(b) = ^(6)(o) = ^ ?&)(«) an< ^ cons fo er tne following decomposition based on the effective metrics 



-(a)(6) - ^(6)(o) - ^(6)( 
5(a) 



C S(b) - r (a)(6)5(c) ■ ( 61 ) 
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(a)(6) 10 "m^Lf^u-^L*^, Wo)(m — C lfc G( a )(f,), and one can choose tensor independent on group 

nik f < ~.ik 



The tensor C^ VM is multiplicative, i.e., C 
index and put 



\a){b) - U (a)(b)9 

This choice can be motivated by the case of parallel fields 



^.tk ^.tk ™ 

9( C ) - 9 9(c) , 



for which the following representation is possible 



T {a)(b)1(c) - G {a)(b) ■ 

The most general decomposition can be much more sophisticated, and we do not consider it here. 



(62) 
(63) 
(64) 



B. Effective metrics for the five-parameter EYMH model 

1. General case 



In the framework of the five-parameter non-minimal EYMH model the macroscopic tensor has multiplicative struc- 
ture, i.e., C^™^) = < ^ lkmn G( a )(b)i the color metrics g 1 ^}^ can be chosen independent on the group indices. One can 
say that in this case there is no "group multi-refringence" , but the spacetime birefringence exists. Then, following 
[46| we represent the tensor Cv^wS (|18p by the decomposition 



r^ikmn / < f 

C (a)(b) ~ b W(6)^ { 



g im(A)gkn(A) _ in(A) g km(A) 



if/ 



i(A) 



jim(B)\/gkn(A) _ gkn(B)\ _ fgin(A) 



in(B) \ fgkm(A) _ n km(B) 
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>]} 



by choosing 



G(A)(A) + G{B)[A) = 1 ) G{A){B) + G(B)(B) = 0> 7 = G^)(B) = G{B)[A) 



(65) 



(66) 



Here the indices (A) and (B) are no longer arbitrary: they are fixed in order to simplify the decomposition of the 
constitutive tensor [46|. We will refer below to this choice as to A-metric, g lm ( A \ and B-metric, g un ( B \ as well as 
to the corresponding A-wave and B-wave. It is reasonable since generally the A-wave and B-wave are not obligatory 
"ordinary" and "extraordinary" , as in classical optics. 
With this decomposition the A-tensor, g tk ( A ) , can be readily written as (see [46| for details) 



ik(A) = jjirjk + ^ £ ik 



The B-tensor, g jfc ( B ), can be represented in terms of eigenvectors X^yX'Ly and X^ of the permittivity tensor e 
PH - These eigenvectors form the triad and possess the properties 



(67) 

k 



9ikX\ a )X(p) = - <5 (q)(/3): 9ikX\ a )U - , 



(68) 



A (1) A (1) + A (2) A (2) + A (3) A (3) - ~9 +UU — — ZA 



7 k 



(69) 



e klm U A (a) A (/3) — e (a)(/3)( 7 )^t7)' 



(70) 



where (a), (/?).... = (1), (2), (3) are tetrad indices; there is a summation over repeating indices. The permittivity 
tensor can be decomposed according to 



£ tk _ -2je (a )X ( * a) -X(* a) , £fc -£(i) +£(2) + e (3), 

(a) 



(71) 
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where the terms £( a ) denote the eigenvalues, corresponding to the eigenvector XJy Since the tensor e lk is orthogonal 

to the four- velocity vector XJ % , the corresponding eigenvalue £( ) is equal to zero, and the velocity does not appear in 
this decomposition. 

The structure of the B-tensor, g im i B \ is much more sophisticated in the general case. Keeping in mind an application 
of this model to cosmology, we present here only one example, when the tetrad components of the impermeability 
tensor (/•* )(a)(/8) are diagonal, i.e., 

(M _1 )(a)(/3) = - — S (»)(f3) > ( 72 ) 
and the tensor of cross-effects vanishes. In this case the B-tensor (see j4fj) has the form 

I (3) 

irn(B) Am{A)^ tl / M(i)M( 2 )M( 3 ) ^ £ {a) ■ 

V 7 (a)=(l) 

where 

M (a) = l ^ . (74) 

M(a)M£(l)£(2)£(3) 



2. Uniaxial case 



Below we consider the applications of the model to the spacetimes with uniaxial three-dimensional subspaces and 
simplify the structure of the A- and B-tensors accordingly to this case. Let the direction of the privilege axis be Ox 3 
and Mm = Mr 2 \ = 0. Then one obtains 

e± = e\ = e\ = 1 + 2Tl u lt , e 3 = e\\ = 1 + 2U 3t 3t , (75) 



— ^^- 1 )l = (^ 1 )l = l + 2TZ 23 23 , (^ 1 ) 3 ^ — = 1 + 2TZ 12 12 , (76) 
-=S||£±M-L, - = 1- — . (77) 

H 7 £_LM|| 



The A- and B-tcnsors take the form 



fHA) = ulxjk _ J_ (xi X k + X f 2) Xf 2) ) - -±-XinXL , (78) 



- /( v L (l)^(l) ^(2)^(2) J - 7^7^(3)^(3) ' 



9 lk{B) = - ^ ( x U x h + x h x w) 7^I x U x k > ( 79 ) 

providing the relation (|65p . In other words, they give the associated metrics for the non- minimal EYMH model. In 
the WKB-approximation the Yang-Mills equations are satisfied when 

ff ife(A Wfe = 0, or 5 lfc ( B W = 0, (80) 

(see Appendix [Bj. Thus, these associated metrics can be interpreted as the color ones. This means, in particular, that 
the propagation of a test Yang-Mills wave in the non-minimally active spacetime can be described by two dispersion 
relations 



2 _ Pi , P \\ 2 _ V\ Pj\ rfm 

(A) - — — + — — > - — + > ( 81 i 



e±fi± v ' e_L/i|| e±(i± 

where 



oj = Pl U l , p 2 ± = -(pip 1 + p 2 p 2 ) , pjj = -p 3 p 3 . (82) 
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When the wave propagates in the longitudinal direction, i.e., p± = 0, or in the transverse one, i.e., pn = 0, the 
corresponding phase velocities can be defined as 

V|| = -, V ± ^— . (83) 
P\\ P± 

For A-wave and B-wave they are, respectively, 



V (A) UB) 1 _ 1 + 2^3 



Vf) = ^- = 1 / l + 27e >, Vf) = ^- = 1 / l + 27el >. (84) 



These quantities depend on spacetime metric, and its derivatives and convert into one, when the non-minimal coupling 
constants vanish, i.e., qi = q 2 = qz = (let us repeat, that in units used here the speed of light in the standard 
vacuum is equal to one). The color-acoustic metric can be represented in the form analogous to (|78p: 



g ik = A 2 

where 



rrirr/e ^" / -yi -y-k , \ri -y-k \ -y-i -y-k 

U U - { A (1) A (1)+ A (2) A (2)) - ^2 A (3) A (3) 



(85) 




^ = i+^- - ^/tt^' °ii = Vi+^I' (86) 

and cr_L and a\\ can be interpreted as phase velocities of a transversal and longitudinal waves, respectively. 

V. COSMOLOGICAL APPLICATION OF THE EFFECTIVE METRIC FORMALISM 

A. One-axis Bianchi-I model 

Relativistic cosmology needs effective metric paradigm since it uses observational data obtained by analysis of 
the directional distribution and frequency properties of incoming photons. Propagating photons are influenced by 
numerous regular and stochastic phenomena, and we could try to identify this influence, if we reconstruct the effective 
dielectric permittivity and magnetic permeability of the regions, where the photons passed. Effective metric formalism 
gives mathematical grounds for such a reconstruction. 

Let us use the metric 

ds 2 = dt 2 - [a 2 {t) (dx 2 + dy 2 ) + c 2 (t)dz 2 ] , (87) 

which is attributed to the well-known Bianchi-I cosmological model with two equivalent spatial directions, Ox and 
Oy. In [6l[ we presented the exact solutions of the non-minimally extended Einstein-Maxwell model with such metric 
in case when there is a magnetic field directed along Oz. That model can be automatically generalized to the case of 
parallel Yang-Mills field F^' — q^ B{8}5 2 — 5 2 5\) along the line proposed in [fijj. Thus, we can refer to the exact 
solutions, published in [6l|, as to the ones with non-minimally coupled parallel Yang- Mills field of the magnetic type. 
This means that we can use both terms: optical metric and color metric referring to that solutions. The tetrad vectors 
for this metric are 

-y-k ?k -yk rfe ^" -y-k sk ^ yk ?k (qq\ 

X {0) -d , X^-d, — , X (2) -6 2 — y X (3) -d 3 —. (88) 

Based on the symmetry of the spacetime we easily obtain that the non-minimal susceptibility tensors 72- mn and 5ft m 
correspond to the case of uniaxial symmetry. The non-vanishing components of these tensors are the following 

^ 1 *i* = ^(49i+3<Z2+2<z 3 )- + i(2g 1 + g2 )- + i(2 (?1 + (Z2 )^+i(4 (Zl + g2 )-- , (89) 
2 a 2 c 2 a z 2 a c 
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-^ 13 i3=^(49i+g2)- + i(2 g i+ g2 )- + i(2 9 i+ ( 7 2 )^+i(4 (Zl +3 (Z2 +2g3)-- , (90) 
2 a I c 2 a z 2 oc 



^ 3t 3t =(2 g i+ g2 )- + ( gi + (Z2 + (?3 )-+ (Z i — + (2gi+ 92 )-- , (91) 



^ 12 12 = (2 9 i+ (Z2 )"+gi^ + ( (Z i+ g 2+g3)^2+(2gi+g 2 )^ , (<>2) 



n 2t 2t = K u lt , ^ 23 23 =ft 13 13 , (93) 

n n \ 

(94) 













■1) 




1-2— 


V a 






ac J 



»i = = (494 + 95) f- + -- ) + 2g 4 - + (2g 4 + g 5 )4 - (95) 
\a a c J c a z 



Kl = 2g 4 (2- + ^)+(2 g4 + (?5 )(- + 2--) . (96) 



a a 2 1 \c a c 



Thus, we are ready to write the effective metrics in terms of permittivity (|T5[) and impermeability (176 



i. Color (optical) metrics 



The A-metric and B-metric are, respectively, 

ik(A) = gigk 1 (gtgk + glgf) (97) 

= gi g H _ ^_ { gi i§ k + gi^ _ -J—**** , (98) 

where em, /in and /ij_ are given by ([75)) and (|75|) with the susceptibility tensors (fB"9")) - (|96p . The relative anisotropy 
of the permittivity/impermeability can be estimated by the following quantities 

£±-e ll ^{q 2 + 2q 3 )[ )+q2-[ , (99 

\c a J a \c a J 

11 / c d\ , d / c d\ . 

=g 2 +(2g 3 + g 2 )- • (100) 

Hi ui \c a J a \c a J 

It is worth noting that q\ does not enter the formulas (|M|) and (|100|) . The isotropization of the Universe leads to the 
vanishing of these differences. 

Special case <? 2 + qs = 

It is clear from (|89p ~ (|93p that in this special case four non- vanishing components of the non- minimal susceptibility 
coincide: 

T^\t = K u u = ^ 23 23 = ^ 13 i 3 , (101) 
and, thus, the permittivity scalars are linked by 

ej_nj_ = 1 , E||//|| = 1 , (102) 
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which provide the relations 



(A) _ _ 



vp = i , v| A) v| s) = 1 , 



II 



The color metrics in this case can also be rewritten as the functions of only one parameter, say, v 2 



ik(A) _ Tjijjk _| 

V 2 



1 I x\ 3) x\ 3) , 
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(103) 

(104) 
(105) 



The structure of these color (optical) metrics shows that the waves propagate in the longitudinal direction with the 
speed of light in standard vacuum. As for the waves propagating in the orthogonal direction, the phase velocity of 
one of them is less than speed of light in vacuum, the second wave being superluminal. We deal with birefringence. 



2. Color-acoustic metric 

Using (|94|) -([96 |) one can write the color-acoustic metric in the form 

1 



where 



f k = A 2 {t){5l5l- 



a?{t)a\{t) 



1 pxk \ 



























-2q 4 




h2— 










a c 



(106) 



(107) 



al(t) 



The eikonal equation for the scalar fields reads 



Here the following definitions are introduced: 



1 - (2?4 


+ *>)(2f- 


h f ) - 2 Qi 




h 2^ 

a c 




l-2 94 


2 ! + (ff 


- (2qr 4 + 95) (| - 


a cy 




1 - (2<?4 


+ 96)(2f- 


f §)-2<Z4 




a c 




1 - (4g 4 + 


9s)(f + |- 


-)-294f- 


- (2g 4 - 


K«0( 


A) 2 

a/ 



~g lk P l P k = m 2 



£ 2 = m 2 A- 2 + fjfo|y 2 + P|crJ 2 , 



£ 2 = P t P l , P\ = -{PiP 1 + P 2 P 2 ) , P,f = -P 3 P 3 , 



(108) 

(109) 

(110) 

(111) 
(112) 



where Pj is particle four- momentum, £ is its energy. When a scalar particle moves along longitudinal direction, i.e., 
Pi = 0, its three- velocity is 



»||(t) = = <7||Vl-mM- 2 f- 2 , 



(113) 



thus, for high-energy scalar particle {£ ^> m) the quantity a\\ gives asymptotic longitudinal velocity. When P|| = 0, 
one obtains, respectively, 



V±(t) = ^r = <J^l-m 2 A- 2 £- 2 , 



(114) 
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thus, <7_l is some asymptotic transversal velocity. The ratio 



^|| ft) = CT||(0 
V±(t) (T±(t) 

does not depend on the particle energy £ and is predetermined by the values of the functions a(t) and c(t) and their 
derivatives, as well as by the values of the non- minimal coupling parameters 94 and 95. Note, that when 95 = and 
the derivative coupling is absent, this ratio is equal to one. 



1- (494 + 95) (| 



a c 

a c . 



- 2q4 - {2q A + q 5 ) (§) 



l-2« 4 



2| + (t) -(294 + 95) 



222 

a c) 



(115) 



Exactly integrable example 

Let us extract from [6l| one of the exact solutions of the non-minimally extended Bianchi-I model with magnetic 

a 2 ( t ) 

field ($ = 0). This solution is characterized by the following features. The magnetic field is B(t) = B(to) g a(°v • It 

can be replaced in our model by the "parallel" Yang-Mills field of the magnetic type with = B(to)a 2 (to)5^y 
The longitudinal V\\ and transverse V± pressures of the matter are connected with the energy density W as follows: 
V\\ = —V\_ = —W. The solution for the a(t) is of de Sitter type 

a(t) = a(t )e H(t - to) , A = 3H 2 , (116) 

with the constant H given by 

2 2W(t ) + B 2 (t ) 

H = 1094i? 2 (io) ' (117) 

the coupling constants being linked by the relation 691 + 492 + 93 = 0. The solution for c(t) is 

C (i) = C (t )e H( *-* o)i tz ; i -, a= Kqi B 2 (t ). (118) 

(1 - a) 

Clearly, this non-minimal cosmological model is non-singular for t > to, if a < 1, i.e., the first coupling parameter 91 
satisfies inequality < 91 < ^gjrc^ ■ For such model the relative anisotropy can be characterized by 

e e Mq 2 +Aq 3 )[2W{h) + B 2 {t )} 

£± " £|l = 5 [a - e**(«-«o)] ' (119) 



J_ _ J_ = 2 K (q 2 -2q 3 )[2W(t ) + B 2 (t )} 

M|| H± 5 [a - e 4ff (*-*°)] ' 1 J 

The anisotropy of the permittivity disappears exponentially at < — > 00, i.e., when Bianchi-I model isotropizcs. Mention 

that e± = £11, when q 2 + 4g 3 = 0, as well as, — = — , when q 2 = 2q 3 , nevertheless, it can not occur simultaneously, 

11 - /1 1 1 [j,^ 

since when q 2 = 93 = 0, 91 also vanishes and the non-minimal model degenerates. Analogously, the color- acoustic 
anisotropy is characterized by the relation 

1 - ae-^(*-*o) - H 2 [3(494 + 95) + (95 - I2 qi )a e - 4 H ^] 



a±J [1 - 3H 2 {4q 4 + 95)] [l - ae-^C*-*")] 

Of course, at t — > 00, this quantity tends to one, ^i- — > 1, as it should be. 

The analysis of the expressions (|84|l for the given model shows, that for generic q 2 and 93 there are, in principle, two 
critical moments of time, when the phase velocities vanish and the A-wave or/and B-wave stop. At the first moment, 
t\ , the phase velocities of both waves propagating in the longitudinal direction, as well as of the A-wave moving in 
the orthogonal direction, vanish. It is possible, when 



^ - >■ ■■■;„>-{ ,: -""'x i ""i 



18 



if the argument of logarithm is more than one. At the second critical moment, t%, the B-wave propagating in the 
transverse directions has vanishing phase velocity, it is possible, when 



if the argument of this logarithm is also more than one. In its turn, we can choose the parameters a, q 2 and (73 so that 
both optical metrics are regular at each time. Let us illustrate this possibility by the example, when q 2 = — 2qi < and 
<73 = 2</i > 0, satisfying the basic condition 6q± + 4^2 + 93 = 0. For this model the waves moving in the longitudinal 
direction have the phase velocity equal to one, as for the waves propagating in the transverse direction, they are 
characterized by the following reduced formulas 

V (A)\ 2 _ (1 - 4 9l g 2 ) - o(l + 4 gl i? 2 )e- 4 ^-M V (B) V (A) { , 

x ) (1 - A qi m) - a(l - 20giH 2 )e-«'(*- f o) ' x x y ' 

Taking (|117|) into account, we obtain that the square of the phase velocity (1124)) is positive for t > to, when 

2 3 W{t Q ) (5a -2 3 -7a] , . 

5 <a< 7' ^M <mm {^'TT^|- (125) 

Analogously, the longitudinal and transversal color-acoustic waves stop, when, respectively, 

*s = *o + tef ln — ou2^„ : - \ r ' ( 126 ) 



4ff I l-3H 2 (4q 4 + q 5 ) 



t;=«0 + ^7lna. (127) 

4:11 



The color- acoustic metric degenerates when A(t^) = 0, where 

l_ f ojl -ff 2 (12 g4 + ll g 5 )] 
4ff n \ 1 - 3ff 2 (4 g4 + <?s) 



Let us mention that the color-acoustic metrics does not exist initially, if the non-minimal parameters are coupled by 

10gi = (4<? 4 + g 6 )[2W(<o) + B 2 (t )} . (129) 

To avoid the singularities in the color-acoustic metric one can put, e.g., q$ = — 4 g4 < 0, providing the expressions 
(fTTJ7|) and (fT2T|) to be positive for t > to- 

B. Spatially flat FLRW model 

FLRW model with k = (i.e., spatially flat model) can be obtained from the previous model, if we put c(t) = ait). 
The formulas for the A- and B- metrics can be converted into 



gik (A) = g ik(B) = § i 5 k _ J_ (fak + 5 i 2§ k + ftgk) t (130) 

where 



Sfia 2 



: = I — 2 (>!(/i - 2ii 2 q 3 )--2(3 qi + q 2 )(-) - (131) 

a V a 



f i=l-2('3q 1 +q 2 )--2(3q 1 + 2q 2 + q 3 )(-) . (132) 



a \a 

Effective refractive index n.( e g) is given by 



2 2 _ 2 W ~ 2 ( 3 gi + 2 g2 + qsYda - 2(3 gi + g 2 )a 2 ] 

n (eff) _ e^a -a ^ _ ^ + ^ _ ^ + ^ + ^ . 2] . (133) 
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The effective refractive index coincides with the standard one, i.e., ri( e g) = a(t), first, when q 2 +(73 = 0, second, when 
j = Hq = const. When the law of isotropic expansion is the power-like one, say, a(t) = a(to) (^j^j > then (|133[) 
reduces to 

n \eS) _t 2 -2uj [w(6gi + 3g 2 + q 3 ) - (3g x + 2 92 + g 3 )] 



a 2 (t) t 2 - 2lu [w(6gi + 3g 2 + ©) - (3gi + ga)] 

At the moment t* 



(134) 



f = [w(6gi + 3(Z 2 + q 3 ) - (% + g 2 )] (135) 

the color (optical) wave stops, since the phase velocity Vp/, = ° vanishes, of course, the expression in the square 

root has to be positive. There is a number of sub-models in which the singularities do not appear. For instance, when 
qi = —q, q 2 = 2q, q 3 = —q, q > (this model has been considered in [23|, [25| for the spherically symmetric static 
case) and u > 1, the expression (|134l) is positive for each time. Color-acoustic metric reads 



~ik 



A 2 (t) 



(136) 



A 2 (t) = l-3(2q 4 + q 5 )-~6q 4 (-) , (137) 
a \a J 

« 2 (t\ n*(t\ „ 2 (t\ l-3(2 g4 + g 5 )|- 6 g4 (g) 2 

l-(6<Z4 + e?5)§-(6<Z4 + 2g 5 )(g) 

The quantity cr differs from one if and only if q§ ^ 0, i.e., the derivative coupling of the Higgs fields with curvature is 
present. Finally, this color-acoustic metric is non-singular, if, for instance, 55 = — 2q 4 > and u> > |. 

C. De Sitter model 

We consider now the model with positive curvature, K > 0, and reduce the metric to the de Sitter form [62j 

ds 2 = dt 2 - exp{2VKt}(S Q pdx a dx! 3 ) . (139) 
For this spacetime the Riemann, Ricci tensors and Ricci scalar, respectively, take the form 

Rikmn = ~K {gim9kn ~ 9in9krn) , Rtk = SK g ik , R = -12K . (140) 

The tensors Hikmn (fl6]) and Uik (fT7|) yield 

TZ-ikmn = -K(6qi + 3q 2 + q 3 ) (g im gkn - gin9km) , 3?* = -3-ftT(4<7 4 + q 5 )g ik • (141) 

Let us mention that, when we deal with de Sitter model, only two combinations of the coupling parameters 
Ai=6qi+3g2+<?3 an d -^4=94+2-^5 (see (|2~4"|) ) enter the susceptibility tensors IZikmn and 3?^. In this subsection we 
use the parameters Ai and A4 instead of qi,...q5 in order to shorten the formulas. Then, the H ifc tensor and the \& m 
vector simplify significantly, and the equations (|28[) and (|30l) convert, respectively, into 

(1 - 2X 1 K) D k F lk = G 2 (1- 12X A K) $] , (142) 

(1 - l2\ A K)b k D k ® = -m 2 <f>. (143) 
Thus, when Ai < and A4 < the color and color-acoustic metrics are simply proportional to the spacetime one 



g (A)ik = g (B)ik = g »k = ^/i_2\ 1 K g lk , (144) 
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~(A)ik _ ~(B) lk = _ UXiR g ik ( W5 ^ 

One can redefine the coupling constant Q and the mass m 

G ^ G * = g \fi^S-> m ^ m * = m frS^K> (146) 

so that the Yang-Mills and Higgs equations take the standard (minimal) form, but the constants are non-minimally re- 
defined. In the WKB-approximation the conditions g* lk PiPk = and g lk piPk = are in fact equivalent, thus, the mass- 
less test particles propagate in the "non-minimally active" de Sitter model as in standard de Sitter one. This fact has 
been emphasized above during the analysis of the FLRW model. The special case, when 2K (6qi+3q2+q3)^l—2XiK 
and 3K (4q4+q^)=l=12\4K , related to the so-called hidden Yang-Mills and Higgs fields, is discussed in [211 ]. 



VI. CONCLUSIONS 



1. We developed the concept of effective metrics and introduced the associated, color and color-acoustic metrics, 
attributed to the Einstein- Yang-Mills-Higgs model. The procedure of reconstruction of the effective metrics consists 
of two steps. The first one, the reconstruction of the associated metrics, is based on an original multi-metric repre- 
sentations (|54p and (|61[) of the constitutive tensors for the Yang-Mills and Higgs fields, respectively. The second step, 
the identification of the associated metrics with color and color-acoustic ones, is provided by the WKB-analysis of the 
master equations and by the well-known analogy with optical metrics. 

2. Curvature interactions between gravitational, gauge and scalar fields are described in terms of five-parameter 
non-minimal EYMH model, and the basic ideas concerning associated, color and color-acoustic metrics are realized. 
The color (see ([78)1 . (|79[) with ([75)1 . ([75)1 ) and the color- acoustic ([55)1 metrics for the uniaxial models are reconstructed 
explicitly. 

3. The formalism of effective metrics is applied to Bianchi-I, FLRW and de Sitter cosmological models, and the asso- 
ciated, color and color-acoustic metrics are represented explicitly in terms of metric coefficients of the corresponding 
spacetimes and their derivatives. Phase velocities of the corresponding non-minimally coupled color (|84[) and color- 
acoustic (|86[) waves are the main distinguishing features in these expressions. These phase velocities and their ratios 
are calculated for all models (see, e.g., (I108p . (|109[) . (|115p . (|121|) . (|124|) 1. Thus, we obtained the tool for a qualitative 
analysis of the problem of propagation of test color, optical and scalar waves in the "non-minimally active" vacuum, 
interacting with curvature. 

4. The analysis of the presented models allows us to conclude that for generic set of coupling constants q±, q%, . . . q$ 
the color and color- acoustic metrics are singular. This means that generally the phase velocities of the corresponding 
waves can be equal to zero or infinity at specific time moments (see, e.g., (j!22p . (|123p . (I126p . (|135p h The first case 
relates to the stopping of the corresponding wave. The second one is associated with the appearance of the analogs 
of trapped surface s, appa rent horizons or event horizons, which are discussed in detail in the literature devoted to 
analog gravity [H, [37ll38j . 

5. In its turn we present the examples of the sets of coupling constants q±, q2, ■ . . 95, for which there are no singularities 
in color and color-acoustic metrics. Such (regular) models predict, in particular, that the phase velocities of the waves 
essentially depend on the time moment, on the direction of propagation and are very sensitive to the choice of the 
values of the free parameters of the model. These conclusions are in agreement with the ones made in numerous 
papers, devoted to the variation of the speed of light under the influence of curvature interactions in electrodynamics 
(see, e.g., [g3| - [(39| and references therein). These conclusions require to attract a special attention to the analysis of 
the birefringence and time delay effects in observational astronomy, since the photons registered simultaneously are 
not necessarily emitted at the same time. We hope to devote a special paper to the study of such problems. 



APPENDIX A: SYMMETRY OF THE CONSTITUTIVE TENSORS 

1. Tensor 

Model under discussion (see (fl~5)l and ([To]) ) is characterized by the tensor symmetric with respect to indices 

(a) and (b), as well as symmetric with respect to transposition of the pairs of indices ik and mn. In the general 
case this symmetry can be broken, and the tensor Cy*™" can loose the symmetry properties CJ k V}^\ = C % $?£. and 
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^(a){b) = ^{a){b) separately. For £7(1) gauge invariant electrodynamics the antisymmetric part of C lkmn : 

/-yikmn ( /~iikmn /~imnik\ ( a i \ 

^(skewon) = ~° 

introduced by Hehl and Obukhov (42j and indicated as skewon, disappears from the Lagrangian. In the SU(n) 
symmetric Yang-Mills-Higgs theory (n > 1) the skewonic terms can appear in the Lagrangian explicitly, since the 

strength field tensor has a group index with two and more values, and indices (a) and (b) can be in the 

antisymmetric combination. Generally, we have the following decomposition of the second term in the Lagrangian 

/-likmn Z7'( a ) Tp(b) ^ r\ikrnn 77i( a ) jp(b) _±_ ^ ^ikrnn T?{ a ) 17(b) ( a o\ 

where the tensors and S^™) possess the following symmetry of indices 

r\ikmn r\mnik r\ikrnn "ifcmn '-j ninth _ ^mnik / a q\ 

U (a)(6) - U (a)(6) - U (6)(a) > -(a)(6) ~ -"(a)(6) ~ "(6)(a) ■ 

The term explicitly represents the Yang-Mills skewon, if we use the terminology introduced by Hchl and 

Obukhov [42]. In this paper we do not consider such a model and assume that = 0. 

The symmetric part of the tensor of linear response S|„)?R can be decomposed using the velocity four- vector U % in 
the same way as in electrodynamics: 

eikmn / 'im jjkjjri in i-j-krj-m , kn rjijjm ~km TT^TT n \ 
(a)(6) — 2 \ £ {a)(b) U U ~ £ (a)(b) U U + £ (a)(b) UU ~ £ (a)(b) UU J 



ik ( —l\la ^mn , ik ijjm in jjn im \ . rnnljn Ik 

V ill" i(a)(b) 7 7 s +? 7 il^ ^(a)(6) U ^o] 1(6) I + Vl ( U V {a 



Ik TT^.M \ 

(a)(6) - U V {a){b)) 



(A4) 



Here £f!Xiu\ and (/i /('vm are ana logs of the tensors of dielectric permittivity and magnetic impermeability, respec- 



-(a)(6) 



tively, and K a \( b \ is an analog of the tensor of magneto-electric coefficients. These quantities are defined as 



■■( a )(b) 



= 29 



ikmn 
(a)(6) 



u k u n 



/(a)(6) - 2 ' ik ( a )( b ) lmn ' 



(A5) 



"(a)(6) - ^ zfc U (a)(6) U n ~ L/ fc U (a)(6) ? 7in ■ l At) J 

Tensors 77 m „i and rf kl are the antisymmetric tensors orthogonal to XJ % defined as 

Vmm = e mn isU s , if 1 ee e lkls U s . (A7) 
They are connected by the useful identities 

-V ikp V m n P = §%lUiU s =A i m A k n -A i n A k m , \^f kl r\kim = -& ms UiU s = -A l m , (A8) 

where <^i s and <5^ s are the generalized Kronecker tensors [Hj]. The symmetric projection tensor A lk is defined as 
A lk — g lk — U l U k . The tensors £ („)(m an d (i" -1 )^)^) are symmetric, but v u\(b) * s generally non-symmetric with 
respect to the spacetime indices. These three tensors are orthogonal to U\ 

ef a)(b) U k =0, (v-% m U k = 0, ^ )(b) ^=0 = ^ )(b) C4, (A9) 

and possess 21 independent components for any fixed set of the indices (a) and (b). 

To complete the analogy, we introduce two four- vectors of excitation DVx and ^\ a y as well as, two four- vectors of 

the field strength £ ) > and B\ > by the definitions 

o (a) (a) J 

-DU=H%U k , H\ a) =H*$U k , £W=F$U k , B^=F* h {a) U k , (A10) 
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where the four- vector of the velocity U l is normalized by unity: U l Ui = 1. These vectors are orthogonal to the velocity 
four- vector U l : 

V\ a) U i = Q = £ < f ) U\ W^Ui = = B\ a) U\ (All) 

and form the basis for the decomposition of the Fmn and TiTff tensors: 

Fmn = £$U n - £ ( n a) U m - e m tU s B\ a) , K$ = V^U n - Vf a) U m - e mn ls U s H\ a) . (A12) 
When we deal with the five-parameter EYMH model, the direct calculations using (|26p show that 

e$ (5) = [A im + m ikmn U k U n ) G [a)(b) , (A13) 

= G W(*>) i APq - 2 *^* plqs UiUs] , (A14) 

"WW = -K^U^G^ . (A15) 

Thus, one can state that the non-minimal coupling of the gravitational and Yang-Mills-Higgs field effectively changes 
the properties of the vacuum linear response. Particularly, the tensor of non-minimal susceptibility TZ^ determines 
the variation of color electric- type permittivity, the corresponding double dual tensor *JZ* plqs is responsible for the 
changes in the color magnetic-type impermeability tensor, and non-vanishing *JZ plqs produces the so-called color 
magneto-electric-type properties of the non-minimal vacuum. 

2. Tensor C[ k a)(b) 
The third term in the Lagrangian {T]) also contains two parts 

- c\ k am b^b k ^ = -n% m Di$WD k $<® - ^ K6) A$ (a) M (6) , (Ai6) 

where tensors f2f^ and "J/^w^ are symmetric and antisymmetric, respectively, i.e., 

(«)(*>) ~~ ~~ (&)(«) ' W («)(6) ~ W (a)(6) ~ W (6)(a) • l Ai 

The decomposition of the symmetric tensor nas ^he same structure as a symmetric stress-energy tensor 

fi (a)(b) = A(a)(b)U l U k +J{ a ) {b )U k +J k a ) {b) U l +S( )(6) , (A18) 

where JL^tfy and BiK,^ are orthogonal to the velocity four- vector U k . For the non-minimal vacuum we obtain 

-4(a)(6) = G {a){b) (1 + q 4 R + q 5 R ik UiUk) , (Af 9) 

J(a)( b) = G (a)(b) q 5 (R im - WR mn U n ) U m , (A20) 

B$W = G(a)( b) [(1 + qiRW' + q 5 R mn Ka^ k r] ■ (A21) 

The skew-symmetric part ^(a)( ) nas the decomposition analogous to the Maxwell tensor 

*W(6) = ^(a)(6)^* " M\ am TP ~ V iM Ma m , (A22) 

with 

M \a)(b) = ~- M \b)(a) > A *(a)(°) = -Aj(6)(o) > (A23) 

and introduces the scalar skewons by the analogy with the terminology of Hehl and Obukhov. 
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3. Tensor C l]kl mnpq 

This tensor can also be decomposed into irreducible parts using the same procedure as in case of ■> nevertheless, 

we do not present this procedure here. 

APPENDIX B: WKB- APPROXIMATION FOR THE GAUGE FIELD IN THE UNIAXIAL CASE 

In the WKB-approximation the gauge potentials Ar and the field strengths F^f' can be extrapolated as follows 



Pk A\ a) - P iA { : ] 



Pk = v fc * 



(Bl) 



Let us mention that the nonlinear terms in ([5]) give the values of the next order in WKB-approximation, thus, such a 
model of gauge field is effectively Abelian. In the leading order approximation the Yang-Mills equations are reduced 
to 



Cfam Pk Pm =0. (B2) 
Substitution of CffyjX from (J65|) with g lk( - A ^> from $8$ and g lk{B ^ from ([79| with jj, and 7 given by J77|) yields 

9 im(A) P m \g kn(A) P k A^] ~g m ^A^ \g km{A) p k p m = 



1 - » ) PkPm A^ 



(1) A (2) ~~ A 



(2) A (1)J ^ A (1) A (2 



2) - A (2) A (1) 



(B3) 



Projection of this equations onto the velocity four-vector {/, gives the scalar ratio 



(U k Pk ) 



g mn ^p m A^] = [u n A { : ] ) 



km(A) 
9 PkPn 



(B4) 



which is satisfied if, for instance, we use the Landau gauge U n A\f > — and the condition of orthogonality of the wave 
four- vector and amplitude four- vector in the first associated metric, i.e., g^PkA^ = 0. Projections onto the axes, 
given by Xm> an d XJ 3 \ yield, respectively, 



A {a) 



g km(A) PkP m +p\ 2) 



1 



1 



£||/i_L £j_/i|| 



A 



(a) 
(2) 



P{X)P{2) 



1 



|M-L £ J-M|| 



0. 



A {a) 



P(1)P(2) 



1 



1 



£||/X_|_ £j_/X|| 



-.4 



(«) 
(2) 



km(A) 



1 



£||M_L e_LM|| 



= 0. 



.4 



(«) 

(3) 



0, 



(B5) 



where A^l = X k ^A k a \ p(i) = X k ^p k , etc. Taking into account the relation 



g km{B) Pk P , 



g km(A) PkPn 



1 



(B6) 



one can conclude that nontrivial solution of (|B5[) exists, when 



g km(A) PkPn 



g lj(B) PiPj 



= 0. 



(B7) 



Thus, the associated metrics g km ^ and g fcm ( s ) are the color ones. 
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APPENDIX C: ALTERNATIVE DESCRIPTION OF THE PARTICLE MOTION IN TERMS OF TIDAL 

FORCE 

Particle dynamics can be alternatively described using the equation of motion with effective force. Instead of 
equation of null geodesies in the effective spacetime with metric g lk ( a \ where a — A, B, one can write the equation 

n^ = F^\ (ci) 



dr 2 dr dr 



where 



( a ) dx k dx l 
d ~d~7~di 



F^nJ-JH-f?-, nfEn r rf. (C2) 



IT, and r fei are the Christoffel symbols for the real and effective spacetimes, respectively. The quantity n fc; , the 
difference of the Christoffel symbols symmetric with respect to indices k and I, is known to be a tensor, thus, the 
quantity F 1 ^ is a vector. Since we consider the interaction of Yang-Mills and Higgs field with spacetime curvature, 
we can indicate this force as a tidal one. The tidal force F 1 ^ is quadratic in the particle velocity four-vector and is 
predetermined by the structure of the tensor nlj . For Bianchi-I model this tensor has the following non- vanishing 
components: 

n n = n 22 = 2 J t i a l 1 "^^)] . n n = n 22 - g J t i a i 1 -^^)! . 



where the quantities en, e_L, /Un and /Ltj_ are defined in (|75p and (fTS)) with the susceptibility tensor components from 
(155)1 - (| 9 3 1) . Thus, the force F l(Q ) consists of derivatives of a(t) and c(i) up to the third order. 
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